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NOTATION 


i 


A  Wave  single  amplitude 

Cp  pressure  coefficient,  (p-po)/  l/2pV,2 
Cpmin  Minimum  pressure  coefficient 

D  Propeller  diameter 

H  Significant  Wave  height 

K»  Thrust  coefficient,  Kt  =  T/(p  n2  D4) 

J»  Advance  coefficient,  V,/(nD) 

r  Local  propeller  radius 

R  Ship  resistance 

T  Propeller  thrust 

4>i  First  harmonic,  time  domain,  phase  angle  at  the  propeller  plane  due  to  seaway 
6  Angular  position  at  propeller  plane,  degrees,  0  @TDC 
Vt  Tangential  velocity  of  the  propeller  wake 

Vx  Axial  velocty  in  the  propeller  plane 

Vr  Radial  velocity  in  the  propeller  plane 

Vs  Ship  Speed 

Va  Volume  average  effective  axial  velocity  at  propeller  plane 

Subscripts 
ave 

c 

f 

peak,  ave 
cw 
min 
max 
t 
x 
r 

SS 
1 


Circumferential  average 
captive,  waves,  but  no  hull  motion 
free,  waves  and  hull  motion 
Average  peak 
Calm  Water 

Value  producing  minimum  loading 

Value  producing  maximum  loading 

Tangential 

Axial 

Radial 

Sea  State 

First  spatial  harmonic  at  the  propeller  plane 
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ABSTRACT 

A  computational  procedure  has  been  developed  to  predict  the  influence  of  ship 
motions  in  waves  on  propeller  cavitation  inception.  A  test  case  has  been 
performed  on  the  DD963  hull  form  and  propeller  predicting  the  influence  of 
various  sea  states  on  propeller  cavitation  inception.  Simple  modifications  of  the 
propeller  geometry  were  made  to  assess  the  potential  of  improving  propeller 
cavitation  performance  with  consideration  of  design  trade-offs. 

ADMINISTRATIVE  INFORMATION 

This  report  is  submitted  in  partial  fulfillment  of  milestone  2,  Propulsor  Design  Methodology, 
of  task  4,  Effects  of  Sea  State,  in  the  Advanced  Propulsion  Systems  Project  (RB2 133-MS  1),  in  the 
Surface  Ship  Technology  Program  Plan  (Program  Element  0602121)  for  fiscal  year  1995.  The 
work  described  herein  was  sponsored  by  the  Office  of  Naval  Research  (ONR  334)  and  performed 
by  the  Carderock  Division,  Naval  Surface  Warfare  Center,  Code540  under  Work  Unit  Number  1- 
5060-554. 


INTRODUCTION 

In  recent  years,  there  has  been  increasing  interest  in  the  consideration  of  propeller 
performance  when  operating  in  a  typical  seaway.  In  some  cases.  Navy  propeller  design 
requirements  determined  by  NAVSEA  have  included  propeller  cavitation  inception  with  ship 
operation  in  a  seaway.  An  example  was  in  the  design  of  an  Oceanographic  Research  Ship,  T- 
AG(X)  1(1989)  where  operation  in  extreme  sea  states  was  considered  by  applying  margins  to  the 
still  water  cavitation  predictions. 

With  present  predictive  tools  for  the  motion  of  a  ship  in  a  seaway,  it  was  appropriate  to  develop 
a  numerical  procedure  to  calculate  the  time  varying  flow  field  into  the  propeller  disk  due  to  a  ship 
operating  in  regular  waves  and  calculate  the  effect  of  the  unsteady  flow  field  on  propeller 
performance.  This  procedure  is  potentially  an  improvement  over  less  accurate  estimations  of  ship 
motions  and  wave  orbital  velocities.  The  procedure  used  for  predicting  the  ship  motion  in  regular 
waves  was  the  SWAN^(Ship  Wave  Analysis)  code.  This  code  is  potential  based,  and  solves  the 
complete  3D  potential  problem.  The  code  output  was  interpreted  to  calculate  the  time  varying 
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velocity  at  the  propeller  plane  by  Kim  and  Chevalier^.  The  propeller  performance  was  calculated 
in  a  quasi  steady  fashion  using  a  potential  based  panel  code^. 

The  usefulness  of  an  improved  prediction  capability  is  to  provide  the  designer  the  ability  to 
consider  the  influence  of  the  ship  in  a  seaway  in  the  overall  requirements  of  the  design.  A  sample 
design  exercise  was  performed  to  illustrate  design  trade-offs  that  can  be  performed. 

PREDICTION  OF  PROPELLER  INFLOW  DUE  TO  WAVES  AND  SHIP  MOTION 

The  prediction  of  propeller  inflow  due  to  waves  and  ship  motion  was  performed  using  the  MIT 
SWAN^  code.  The  output  of  the  code  was  modified  by  Kim  and  Chevalier^  to  compute  the  time 
varying  velocity  in  the  propeller  disk.  The  entire  code  was  run  on  a  workstation  with  cases  running 
about  15  minutes.  The  calculation  procedure  was  validated  by  comparing  with  measured  propeller 
plane  velocity  obtained  by  Aalbers  and  van  Gent^.  Both  captive  and  free  model  cases  were 
compared,  with  the  captive  case  correlating  well,  while  the  free  model  case  over  predicted  the 
magnitude  of  die  measured  velocities.  This  implies  that  for  the  free  model  case,  which  is  of  interest 
for  the  actual  design  problem,  the  prediction  method  will  produce  a  conservative  estimate  of  the 
effect  of  the  freely  pitching  and  heaving  ship  in  a  seaway. 

Calculations  were  performed  on  the  DD963  hull  form  operating  in  regular  ahead  seas.  The  full 
scale  ship  length  is  520  feet(162M),  which  was  operated  at  20  kts  in  regular  waves  of  length  equal 
to  ship  length,  with  wave  amplitudes  of  0.6,  1,  and  2.5M,  which  represented  significant  wave 
heights  corresponding  to  sea  states  3,4,  and  5.  In  the  wave  amplitude  range  investigated,  it  was 
found  that  the  perturbation  velocities  in  the  propeller  plane  were  linear  in  magnitude  with  wave 
amplitude.  Consequendy,  future  calculations  would  require  only  one  wave  height  calculation,  with 
assumed  linear  behavior  up  to  at  least  a  wave  amplitude/ship  length  equal  to  0.015. 

VELOCITY  VARIATION  AT  THE  PROPELLER  PLANE 

The  velocity  variation  at  the  propeller  plane  was  represented  as  a  first  harmonic  amplitude  with 
wave  frequency  at  specified  radial  and  circumferential  locations  in  the  propeller  plane.  The  axial 
and  tangential  velocity  was  considered  dominant  in  altering  the  propeller  loading.  Table  1  shows 
the  tabulated  perturbation  velocities  for  the  1  meter  wave  case. 

Figure  1  shows  the  calculated  variation  in  the  first  harmonic  amplitude  of  velocity  as  a 
function  of  angular  location  in  the  propeller  disk  at  r/R=0.57  for  the  free  model  case.  The  axial 
velocity,  V*i,  is  relatively  constant  through  the  propeller  plane,  with  some  increase  in  velocity 
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perturbation  at  the  top  of  the  disk,  presumably  due  to  the  influence  of  the  hull.  The  tangential 
velocity  shows  a  large  variation  with  angle,  which  is  primarily  a  geometric  effect.  The  vertical 
velocity  is  the  primary  lateral  flow  variation,  which  is  also  relatively  constant  through  the  disk. 
When  represented  in  the  cylindrical  coordinate  system  of  the  propeller,  as  a  tangential  velocity  it 
shows  a  once  /revolution  variation.  When  shown  as  the  first  harmonic  amplitude,  it  results  in  peak 
values  at  the  sides  of  the  disk,  90  and  270  degrees. 

From  inspections  of  other  radii  along  with  that  shown  in  Fig.  1,  the  following  assumptions  were 
made  to  simplify  the  representation  of  the  velocity  perturbation  due  to  waves  and  ship  motion. 

1.  The  axial  perturbation  velocity  at  each  propeller  radii  is  assumed  constant  , 
represented  as  the  calculated  radial  average,  (Vx.i)«ve 

2.  The  tangential  perturbation  velocity,  is  represented  as  the  average  of  the 
maximum  values  on  each  side  of  the  disk,  (Vu)pe«k.«ve. 

3.  The  phase  between  the  axial  and  tangential  velocity  perturbations,  is  assumed 
to  be  90  degrees,  thus  permitting  independent  assessments  of  their  effects  on 
cavitation  performance. 

The  table  below  shows  the  amplitudes  and  time  domain  phase  angles  of  the  first  harmonic  of 
velocity  due  to  the  waves  and  ship  motion. 


@0.57R 

(Vxj)ave/V  S 

^l.vx 

(degrees) 

(VuJpeafc.ve/Vs 

*t*l,vt 

(degrees) 

A<|> 

Captive  model 

0.039 

-146 

0.0235 

-58 

88 

Free  model 

0.0516 

-155 

0.049 

-85 

70 

For  the  captive  case,  the  axial  and  tangential  velocity  are  very  close  to  90  degrees  out  of  phase.  For 
the  free  model  case,  this  is  less  of  the  case.  To  simplify  the  analysis  ,  the  velocity  components  are 
treated  independently.  This  assumption  could  reduce  the  peak  influence  of  the  waves  for  the  free 
model  case  by 
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1/(1+  sin(70°))  =  0.75 


When  performing  an  analysis  for  a  specific  application,  then  relative  phase  between  the  two 
components  of  velocity  could  be  easily  accounted  for. 

At  each  radial  location  the  effect  of  waves  and  ship  motion  is  represented  by  a  perturbation  of  the 
still  water  wake  by  (V^O^e,  and  (Vt.Opok.  «ve-  Figure  2  shows  the  variation  of  these  quantities  with 
radius  for  the  1  m  wave  case.  Also  shown  is  the  calm  water  first  harmonic  tangential  velocity.  First 
to  be  noted  is  that  the  1  meter  wave  case  does  clearly  present  a  perturbation  to  the  still  water  case, 
since  the  tangential  velocity  is  about  l/3rd  its  value.  The  increase  in  the  tangential  velocity  for  the 
still  water  case  near  the  root  of  the  propeller  is  due  to  flow  augmentation  from  the  upstream  shaft. 
The  calculation  performed  using  the  SWAN  code  neglect  the  shafting  and  show  no  augmentation  at 
the  inner  radn.  Also  note  that  the  free  model  case  produces  larger  velocity  perturbations  that  the 
captive  case. 

ADDED  RESISTANCE 

The  effect  of  waves  and  ship  motion  on  the  time  average  performance  of  the  propeller  is  a  very 
important  component  of  the  propeller  load  augmentation.  This  aspect  unfortunately  is  not 
addressed  numerically  in  this  procedure.  A  simple  empirical  approach  was  used  to  predict  added 
resistance  for  the  combatant  hull  form  investigated.  Data  from  Lloyd5,  shown  in  Figure  3,  was 
used  to  develop  multiplicative  factors  to  be  applied  to  calm  water  resistance  for  effects  of  wind  and 
waves  in  head  seas.  Figure  4  shows  the  variation  in  added  resistance  for  the  three  sea  states 
considered.  The  table  below  shows  the  multiplicative  factors  used  from  Figure  4.  A  very  large 
effect  of  speed  is  seen  in  the  figure,  with  a  large  increase  in  resistance  at  low  speed  at  high  sea 
states. 


SEA  STATE 

10KTS 

15KTS 

20KTS 

25KTS 

30KTS 

32KTS 

3 

1.131 

1.085 

1.047 

1.033 

1.025 

1.023 

4 

1.335 

1.219 

1.12 

1.085 

1.064 

1.06 

5 

3.29 

2.13 

1.52 

1.35 

1.235 

1.216 

Multiplicative  factor,  M,  where,  Rss  =  M  *  water 


For  reference  purposes,  the  following  table  relates  the  various  wave  height  parameters  investigated. 
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reg.  wave 

amplitude, 

A  (m) 

significant  wave 

height,(m) 

H  =2  A/V2 

Assumed 

Sea  State 

NATO  Sea  State  Standard, 

Significant  Wave  Height(m)  for 

assumed  sea  state 

0.62 

0.88 

3 

0.88 

1 

1.4 

4 

1.88 

2.5 

3.54 

5 

3.25 

CALCULATION  OF  MEAN  PROPELLER  OPERATING  CONDITION  IN  WAVES 
From  the  added  resistance,  calm  water  resistance,  and  the  open  water  performance  of  the  propeller, 
the  operating  advance  coefficient,  Ja,  was  calculated  for  use  in  the  propeller  cavitation  predictions. 
The  procedure  is  as  follows.  Provided  from  powering  data,  EHP,  augmented  by  the  sea  state,  and 
(1-t)  are  specified  for  measured  speed.  The  speeds  at  which  calculations  were  made  were  10,  15, 
20,  25,  30,  and  32  knots.  Thrust  per  shaft  for  this  twin  shaft  configuration  is, 

T  (lb.)  /  shaft  =  (M*326*EHPtotal)  /  (2  Vs  (1-t)),  Vs  in  kts 

To  obtain  the  operating  point,  Kt/Ja2  is  calculated  and  plotted  against  the  open  water  Kt  verses  J. 

Kt/Ja2  =  T/pD2V2 

plotting,  Kt=  (Kt/Ja2)  J2 

Figure  5  shows  a  sample  of  this  graphical  process  for  determining  the  propeller  advance 
coefficient,  Ja,  for  various  speeds. 

A  straightforward  calculation  of  cavitation  number,  c,  was  made  through  the  speed  range  of  the 
propeller. 


a  =  (2gH)  /  Vs2,  where  H  is  the  depth  of  submergence  in  calm  water 
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EFFECTS  OF  SEAWAY  ON  CAVITATION  INCEPTION 
The  cavitation  inception  prediction  method  utilizes  wake  predictions  along  with  propeller  panel 
method  calculations  to  predict  the  blade  pressure  distribution  and  thus  the  minimum  blade  surface 
pressure,  and  the  propeller  cavitation  inception  index,  where, 

ai  “  ■^-'pmin)  calculated  w/  panel  method 


WAKE  INPUT 

The  propeller  inflow  wake  used  in  the  calculations  included  an  idealized  calm  water  (cw)  wake, 
and  the  perturbation  velocities  due  to  waves  and  ship  motion  superimposed  upon  the  calm  water 
inflow.  The  calm  water,  circumferential  average  wake  was  idealized  as  follows, 

Vx  =  1.0,  Vt  =  0.0,  Vr  =  0.0,  for  all  values  of  r/R 

The  maximum  and  minimum  loading  conditions  of  the  propeller  blade  due  to  the  calm  water  blade 
inflow  variations  were  idealized  by  considering  the  flow  inclination,  only,  as  follows. 


^tWcwflinc  ~  Vu(r),  Vt(r)  cwynin  ™"  -Vu(r) 


Where  Vt,i(r)  is  obtained  from  a  standard,  straight  ahead,  calm  water  wake  survey6. 

The  wave  and  ship  motion  induced  perturbation  velocities  were  added  to  the  calm  water  wake  in 
a  fashion  to  create  the  four  independent  extremes  in  propeller  loading  as  follows,  for  the  captured 
(c)  and  free  (f)  model  cases. 

1.  maximum  loading  due  to  axial  velocity , 


VxOOmax.ci' —  1.0  -  (Vx,1(r)»ve)c>f  ,  Vt(r)—  Vr  =  0.0 
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2.  minimum  loading  due  to  axial  velocity , 


Vx(rW=1.0  +  (V*,(r).v.)c,f  ,  V,(r)= -Vul(r)cw,  Vr  =  0.0 

3.  maximum  loading  due  to  tangential  velocity, 

Vx(r)=  1.0  ,  Vi(T)mu,cS=  V«(T)m  +  Vuttp^  C)f  ,  Vr  =  0.0 

4.  minimum  loading  due  to  tangential  velocity , 

Vx(f)=  1.0  ,  Vt(r)  min.cj1  =  -V^i(r)cw  -  Vul(r)  pealc.«ve,  c,f  ,  Vr  =  0.0 


For  the  DD-963  hull,  the  components  necessary  to  construct  the  above  propeller  inflow  wakes  are 
as  follows. 


Vx.dO.ve 

Vt,i(r)  pcak,ave 

r/R 

captive 

free 

captive 

free 

calm  water 

0.37 

0.038 

0.0502 

0.0235 

0.049 

0.178 

0.57 

0.039 

0.0516 

0.0235 

0.049 

0.129 

0.8 

0.040 

0.0569 

0.0240 

0.049 

0.128 

1.0 

0.042 

0.0561 

0.0265 

0.049 

0.120 

PROPELLER  PANEL  METHOD 

A  panel  method  is  used  to  calculate  the  time  varying  propeller  blade  pressure  distribution.  The 
code  used,  PSF104,  calculates  the  steady  blade  pressure  distribution,  assuming  uniform  steady 
inflow.  Both  the  hull  and  wave  motion  is  assumed  low  frequency  relative  to  the  propeller  rotation 
rate,  therefore,  the  effects  of  seaway  is  considered  in  a  quasi-steady  fashion.  The  extreme  wake 
inflows  presented  in  the  previous  section  are  calculated  as  if  the  flow  was  circumferentially 
uniform  and  time  independent. 
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A  further  quasi-steady  assumption  has  been  made  to  simplify  the  analysis.  The  primary 
unsteady  flow  for  the  calm  water  case  is  due  to  the  propeller  shaft  inclination,  causing  a  once  per 
revolution  blade  section  angle  of  attack  variation.  The  reduced  frequency  of  the  blade  unsteady 
load  is  low,  permitting  a  further  quasisteady  assumption  in  the  propeller  analysis  of  the  primary 
shaft  inclination  wake.  The  effects  of  the  shaft  and  struts  are  ignored  which  is  consistent  with 
computation  of  the  inflow  effects  of  waves  and  ship  motion. 

These  assumptions,  along  with  the  wake  details  described  in  the  previous  section,  and  the 
propeller  geometry  were  sufficient  to  perform  the  cavitation  inception  analysis.  Because  the  free 
model  case  was  more  representative,  and  resulted  in  larger  effects  on  cavitation,  the  free  model 
case  was  predicted,  the  captive  model  case  was  not. 

The  panel  method  has  been  demonstrated  to  predict  the  cavitation  inception  performance  of 
blade  surface ,  and  root  cavitation7.  The  prediction  of  tip  vortex  was  not  as  accurate,  but  is  thought 
sufficient  to  predict  trends  in  tip  vortex  inception.  The  paneling  used  to  calculated  the  blade 
pressure  distribution  is  shown  in  Fig.  6.  A  30  by  50  grid  is  used  on  the  blade  as  a  reasonable  trade¬ 
off  between  spatial  blade  resolution  and  run  time  (approx.  10  minutes).  Root  and  tip  cavitation 
inception  were  identified  by  the  minimum  pressure  occurring  at  the  nearest  panels  to  the  hub  and 
the  outer  chordwise  row  of  tip  panels.  Blade  surface  cavitation  was  identified  by  minimum 
pressure  calculated  amonst  the  remaining  panels,  generally  occurring  at  the  blade  leading  edge. 
Calculations  were  performed  over  a  range  of  J,  about  the  operating  J’s  calculated  at  the  calm  water 
condition  and  the  included  added  resistance.  A  typical  panel  method  output  summary  is  reproduced 
in  Table  2,  which  also  shows  the  propeller  input  geometry. 

In  the  present  analysis,  there  was  no  consideration  of  variation  in  the  static  pressure  at  propeller 
plane.  The  head  was  assumed  constant  at  the  still  water  propeller  shaft  axis  value,  independent  of 
waves,  hull  motion  and  propeller  angular  location. 

REPRESENTATION  -  CAVITATION  INCEPTION  CURVES 

Cavitation  inception  curves  were  prepared  resembling  a  traditional  cavitation  performance 
diagram.  The  propeller  cavitation  number  for  operating  points  through  the  speed  range  was  plotted 
verses  the  advance  coefficient,  J.  Also  plotted  was  the  cavitation  inception  curves  for  the  various 
types  of  cavitation,  suction  and  pressure  side  leading  edge  (SSLE,  PSLE),  suction  and  pressure 
side  root,  suction  side  and  pressure  side  tip  vortex(SSTV,  PSTV).  For  this  study,  the  predicted 
inception  speeds  were  obtained  graphically. 
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Fig.  7  shows  the  cavitation  inception  prediction  for  the  effects  of  added  resistance  only.  In  the 
figure  only  the  propeller  operating  points  change  with  increasing  wave  height.  Very  significant 
reductions  in  cavitation  inception  speed  occur  for  the  suction  side  inception  with  the  reduction  in 
advance  coefficients  with  increased  sea  state.  Also  shown  is  a  comparison  of  SSTV  prediction 
from  model  test  data  and  predictions  using  the  method  described.  Obviously  for  tip  vortex,  only 
trends  can  be  predicted  with  this  method. 

Fig.  8a  shows  the  cavitation  diagram  for  1  meter  amplitude  waves  with  free  motion, 
representing  sea  state  4.  Only  the  suction  side  cavitation  is  shown  since  pressure  side  cavitation  is 
not  degraded.  It  is  clear  that  the  cavitation  most  effected  is  tip  vortex,  then  blade  surface  and  lastly 
root  cavitation  Fig.  8b  shows  the  similar  diagram  for  the  SS3,  0.6  meter  wave  amplitude  case. 
Both  curves  show  that  the  axial  velocity  perturbation  has  a  larger  effect  than  the  tangential 
velocity. 

Fig.  9  plots  the  cavitation  inception  speeds,  in  knots,  predicted  as  a  function  of  significant  wave 
height.  As  mentioned  earlier,  the  largest  effect  occurs  with  tip  vortex,  the  least  with  root  cavitation. 
Again,  it  should  be  noted,  that  these  predictions  do  not  include  wakes  from  the  shaft  and  struts, 
which  could  result  in  higher  predicted  inception  speeds  than  is  typical  for  model  or  full  scale  data. 

CONSIDERATION  OF  SEAWAY  EFFECTS  IN  PROPELLER  DESIGN 

With  the  procedures  presented  above,  predictions  of  cavitation  inception  in  a  seaway  can  be 
estimated.  Consideration  of  seaway  performance  can  be  included  in  the  overall  design  of  the 
propeller. 

An  exercise  was  conducted  to  perform  a  propeller  design  iteration  investigating  effects  of 
geometry  changes  on  the  cavitation  performance  in  a  seaway.  A  primary  parameter  effecting 
cavitation  inception  is  blade  thickness.  Fig.  10  depicts  a  variant  of  the  original  propeller  thickness 
to  chord  ratio.  Using  the  same  sea  state  4  propeller  inflow  variation  as  previously  discussed,  panel 
calculations  were  performed,  along  with  calm  water  calculations.  Fig.  1 1  shows  the  blade  surface 
leading  edge  cavitation  inception  diagram  for  calm  water  and  sea  state  4,  due  to  axial  inflow 
variations.  The  modified  propeller  shows  only  a  slight  improvement  of  0.2  knots  inception  speed  in 
calm  water ,  with  a  1  knot  improvement  in  sea  state  4. 

A  noticeable  improvement  is  observed  in  inception  speed  ,  but  there  is  a  design  trade-off  of 
increased  risk  of  thrust  breakdown  at  higher  speeds.  Fig.  12  shows  the  calculated  mid-span  blade 
pressure  distribution  at  high  loading  due  to  axial  variations  in  the  inflow  wake.  At  the  blade  leading 
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edge,  a  reduction  in  the  suction  peak  is  seen.  Mid-chord,  a  1  knot  loss  in  back  bubble  cavitation 
inception  is  shown.  Qualitatively,  thrust  breakdown  will  occur  about  1  knot  earlier  with  the  thicker 
blade. 


CONCLUSIONS 

A  prediction  method  has  been  developed  for  propeller  cavitation  inception  in  a  seaway.  The 
method  is  based  on  potential  flow  calculation  of  the  flow  around  a  hull  in  a  seaway  free  to  pitch 
and  heave.  The  propeller  inception  calculations  are  performed  using  a  quasisteady,  potential  based 
panel  method. 

Calculations  were  performed  for  a  surface  ship  transom  stem,  frigate  type  combatant  hull  form 
operating  in  regular,  ahead  seas.  Wave  amplitudes  corresponding  to  sea  states  3,  4  and  5  were 
investigated.  The  following  conclusions  were  reached. 

1.  The  hull  form,  seaway  calculations  and  the  propeller  panel  method  calculations  were  not 
computationally  burdensome,  and  would  be  appropriate  for  use  in  the  propeller  design. 

2.  Of  the  flow  variations  calculated  at  the  propeller  disk,  the  axial  velocity  variations  were 
somewhat  larger  that  the  tangential  velocities,  and  resulted  in  a  larger  effect  on  cavitation 
inception  speed. 

3.  Of  the  three  forms  of  cavitation  predicted,  the  loss  in  inception  speed  due  to  seaway  is  largest  for 
tip  vortex,  then  blade  surface,  and  least  for  root  cavitation. 

4.  The  prediction  of  inception  speed  lost  as  a  function  of  wave  height  could  be  considered  typical 
for  transom  stem  combatant  hull  forms. 

5.  Propeller  modifications  to  improve  seaway  cavitation  performance  must  be  traded  off  against 
adverse  effects  of  cavitation  at  high  speeds,  such  as  thrust  breakdown. 
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Fig  1.  First  Harmonic  Perturbation  Velocity  in  Propeller  Plane  at  r/R=0.57  Due 
to  Free  Motion  in  a  1  Meter  Wave  Amplitude 
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Fig  2.  Perturbation  Velocities  due  to  1  Meter  Wave  Amplitude 
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SHIP  SPEED,  KNOTS 

Fig  4.  Estimated  Increase  in  Resistance  from  Lloyd-*  data 
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Advance  Coefficient,  Ja 


Fig  5.  Determination  of  Propeller  Advance  Coefficient,  J,  in  Sea  State  3 
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Fig  6.  Propeller  and  Hub  paneling  for  DD-963  Propeller 
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Fig  8.  Cavitation  Inception  Prediction  in  Sea  States  3  and  4 
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Fig  9.  Predicted  Cavitation  Inception  Speeds  for  increasing  Significant  Wave  Height 
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Fig  10.  Modification  of  Propeller  Blade  Thickness 


Fig  1 1.  Predicted  Improvement  in  Cavitation  Inception  Speed  with  modified 
Propeller  Geometry 
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-Cp,  PRESSURE  COEFFICIENT 


FRACTION  OF  CHORD,  0.7R 


Fig  12.  Predicted  Blade  pressure  distribution  at  Maximum  Loading  for  Original  and 
Modified  Geometry 
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Table  1  Calculated  Harmonic  Amplitude  of  Velocity  at  the  propeller  plane  due  to  1  meter 
regular  wave  with  ship  motion 

Table  la.  Captive  Model  Case 


Angle,0 

r/R= 

0.37 

r/R= 

0.57 

r/R- 

=0.8 

r/R= 

1.02 

deg. 

Vr,l 

mm 

Vr,i 

mm 

Vr,l 

him 

0. 

,-0.019, 

113. 

-0.018, 

112  . 

-0.017, 

no.  , 

-0.015, 

107  . 

20. 

,-0.024, 

118. 

-0.025, 

118. 

-0.026, 

118.  , 

-0.028, 

117. 

40. 

, -0.027, 

122. 

-0.029, 

122. 

-0.032, 

122 .  , 

-0.036, 

122. 

60. 

,-0.025, 

125. 

-0.027, 

125. 

-0.030, 

126.  , 

-0.033, 

125. 

80. 

, -0.019, 

130. 

-0.020, 

130. 

-0.022, 

130.  , 

-0.024, 

129. 

100. 

,-0.010, 

139. 

-0.010, 

139. 

-0.011, 

137.  , 

-0.011, 

141. 

120. 

,-0.003, 

-148. 

-0.003, 

-130. 

-0.003, 

-127.  , 

-0.003,- 

-116. 

140. 

,-0.009, 

-74. 

-0.011, 

-73. 

-0.011, 

-73.  , 

-0.013, 

-69. 

160. 

,-0.016, 

-66. 

-0.018, 

-65. 

-0.018, 

-64.  , 

-0.019, 

-63. 

180. 

,-0.021, 

-62. 

-0.022, 

-62. 

-0.022, 

-62.  , 

-0.022, 

-61. 

200. 

,-0.023, 

-61. 

-0.022, 

-61. 

-0.022, 

-61.  , 

-0.022, 

-61. 

220. 

,-0.020, 

-58. 

-0.020, 

-59. 

-0.019, 

-60.  , 

-0.018, 

-60. 

240. 

,-0.018, 

-56. 

-0.016, 

-57. 

-0.014, 

-56.  , 

-0.013, 

-59. 

260. 

,-0.013, 

-51. 

-0.011, 

-53. 

-0.009, 

-53.  , 

-0.007, 

-55. 

280. 

,-0.007, 

-42. 

-0.006, 

-40. 

-0.005, 

-37.  , 

-0.003, 

-31. 

300. 

,-0.003, 

15. 

-0.003, 

12. 

-0.003, 

24.  , 

-0.003, 

32. 

320. 

,-0.006, 

87. 

-0.006, 

79. 

-0.005, 

65.  , 

-0.004, 

52. 

340. 

,-0.013, 

106. 

-0.011, 

101. 

-0.009, 

93.  , 

-0.007, 

80. 
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Table  la  (continued) 
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Table  lb.  Free  Model  Case 


r/R= 

0.37 

r/R= 

0.57 

r/R-0.8 

f 

1.02 

deg. 

Vu 

Vx.1 

♦i 

V*,i 

Him 

mam 

<f>l 

0.0 

0.0502  - 

154.2 

0.0513 

-156.2 

0.0503 

-155.7 

0.0543 

-156.3 

20.0 

0.0500  - 

156.2 

0.0516 

-156.4 

0.0535 

-156.7 

0.0555 

-157 . 0 

40 . 0 

0.0498  - 

156.3 

0.0514 

-156.8 

0.0569 

-155.0 

0.0561 

-157.7 

'60.0 

0.0493  - 

156.4 

0.0503 

-156.7 

0.0519 

-157.4 

0.0525 

-158.8 

'80.0 

0.0484  - 

156.4 

0.0490 

-156.7 

0.0512 

-155.8 

0.0506 

-157 . 5 

100.0 

0.0477  - 

155.1 

0.0496 

-155.2 

0.0476 

-156.8 

0.0457 

-157 . 9 

120 . 0 

0.0465  - 

156.2 

0.0458 

-156.3 

0.0454 

-156.5 

0.0450 

-156.3 

140.0 

0.0459  - 

156.0 

0.0445 

-156.2 

0.0437 

-155.3 

0.0428 

-156.4 

160.0 

0.0457  - 

155.8 

0.0451 

-159.2 

0.0428 

-155.9 

0.0416 

-156 . 0 

180 . 0 

0.0452  - 

154.9 

0.0437 

-155.7 

0.0424 

-155.7 

0.0398 

-155 . 1 

200.0  i 

0.0450  - 

156.9 

0.0435 

-156.8 

0.0424 

-157.0 

0.0413 

-155 . 4 

220.0 

0.0458  - 

155.4 

0.0431 

-153.1 

0.0471 

-152.4 

0.0422 

-155 . 1 

240.0 

0.0439  - 

155.1 

0.0454 

-155.3 

0.0449 

-155.4 

0.0436 

-154.9 

260 . 0 

0.0465  - 

155.3 

0.0460 

-155.0 

0.0458 

-155.0 

0.0460 

-153 . 9 

280.0 

0.0475  - 

155.5 

0.0475 

-155.3 

0.0474 

-155.0 

0.0475 

-154 . 9 

300.0 

0.0481  - 

155.4 

0.0484 

-155.2 

0.0489 

-155.1 

0.0494 

-154 . 9 

320.0 

0.0491  - 

154.4 

0.0498 

-154.2 

0.0504 

-155.3 

0.0510 

-155.2 

340 . 0 

0.0481  - 

152.9 

0.0506  ■ 

-155.8 

0.0517 

-155.7 

0.0527 

-155.7 
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Table  lb  (continued) 


Angle, 0 


deg. 


r/R=0.37 


r/R=0.57 


r/R=0.8 


r/R=1.02 


Angle,© 

& 

=0.37 

r/R= 

=0.57 

II 

00 

o 

r/R=1.02 

deg. 

Vu 

[■B 

Vu 

4>i 

Vu 

■mi 

Vt,i 

mm 

0.0 

0.0145  • 

-124.8 

0.0149 

-130.4 

0.0142 

-123.3 

0.0178 

-130.8 

20.0 

0.0284  • 

-105.0 

0.0302 

-106.3 

0.0324 

-107.8 

0.0348 

-109.3 

40.0 

0.0404 

-96.5 

0.0417 

-97.3 

0.0461 

-100.8 

0.0446 

-99.0 

60.0 

0.0471 

-91.2 

0.0476 

-90.9 

0.0476 

-91.4 

0.0480 

-90.4 

80.0 

0-0483- 

-86.2 

0.0479 

-85.4 

0.0497 

-86.5 

O'.  0469 

-83.3 

100.0 

0.0443 

-81.7 

0.7)435 

-80.4 

0.0425 

-79.0 

0.0415 

-76.5 

120.0 

0.0358 

-76.7 

0.0349 

-75.3 

0.0339 

-73.4 

0.0333 

-72.0 

140.0 

0.0239 

-68.7 

0.0233 

-67.2 

0.0237 

-63.7 

0.0221 

-64.9 

160.0 

0.0109 

-46.3 

0.0108 

-44.1 

0.0106 

-44.1 

0.0108 

-42.2 

180.0 

0.0105 

45.2 

0.0092 

51.9 

0.0078 

50.5 

0.0070 

57.0 

200.0 

0.0221 

78.6 

0.0208 

79.8 

0.0195 

80.6 

0.0184 

81.5 

220.0 

0.0355 

86.7 

0.0337 

86.5 

0.0331 

88.4 

0.0315 

89.6 

240.0 

0.0448 

90.1 

0.0439 

90.4 

0.0430 

90.7 

0.0420 

91.5 

260.0 

0.0492 

92.1 

0.0494 

91.7 

0.0492 

91.8 

0.0493 

91.5 

280.0 

0.0479 

93.9 

0.0488 

93.2 

n  .  0499. 

92.6 

0.0508 

91.8 

300.0 

0.0402 

97.7 

0.0416 

96.0 

0.0433 

94.5 

1T.0452- 

93.0 

320.0 

0.0271 

105.1 

0.0282 

103.1 

0.0315 

102.1 

0.0312 

98.3 

340.0 

0.0112 

125.9 

0.0121 

133.9 

0.0123 

132.2 

0.0127 

129.8 

Table  2  Propeller  Panel  Method  Calculation  -  Typical  Summary  Output 


Received  from  MIT  6/1/93 

Local  Drag  Coefficients  added  1/18/94 

Modified  by  Jim  Bailar  on  or  after  6/17/93 

PROP  4661  J-l.l  VT-  VX+B  UNSTEADY  WAKE  JUNE  6,1995 


ADVCO 

CUT  NBLADE  NC 

NR  MTIP  ITHK 

1.1000  1 

.0000 

5  50 

30  0  8 

R/RO 

P/D 

XS/D 

SKEW  C/D  FO/C 

to/c 

0.300 

1.1650 

0.0096 

2.985  0.1780  0.0000 

0.2360 

0.350 

1.2960 

0.0123 

3.481  0.2100  0.0050 

0.1771 

0.450 

1.4800 

0.0186 

4.810  0.2710  0.0209 

0.1070 

0.550 

1.5660 

0.0266 

6.631  0.3270  0.0267 

0.0691 

0.650 

1.5660 

0.0356 

8.978  0.3740  0.0256 

0.0476 

0.750 

1.4980 

0.0447 

11.895  0.4060  0.0209 

0.0360 

0.850 

1.3810 

0.0532 

15.410  0.4090  0 . 0151 

0.0298 

0.900 

1.3060 

0.0567 

17.403  0.3870  0.0122 

0.0284 

0.950 

1.2220 

0.0596 

19.557  0.3260  0.0094 

0.0279 

1.000 

1.1280 

0.0615 

21.876  0.0000  0.0000 

0.0000 

.... 

SINGLE  BLADE  FORCE  BY  PRESSURE  INTEGRATION 

100  100 

M 

RP/RO 

BKFX 

BKMX  BCL 

BCDP 

1 

0.312 

-0.0045 

-0.0029  -0.0025 

-0.0057 

3 

0.35& 

-0.0020 

0.0059  0.0061 

0.0046 

5 

0.405 

-0.0765 

0.0466  0.0482 

0.0122 

7 

0.452 

-0.2102 

0.1134  0.0915 

0.0176 

9 

0.498 

-0.3787 

0.1975  0.1237 

0.0215 

11 

0.545 

-0.5598 

0.2886  0.1424 

0.0235 

13 

0.592 

-0.7352 

0.3762  0.1493 

0.0235 

15 

0.638 

-0.8887 

0.4499  0.1470 

0.0217 

17 

0.685 

-1.0053 

0.5007  0.1380 

0.0187 

19 

0.732 

-1.0710 

0.5209  0.1243 

0.0149 

21 

0.778 

-1.0751 

0.5070  0.1077 

0.0109 

23 

0.825 

-1.0079 

0.4575  0.0894 

0.0073 

25 

0.872 

-0.8613 

0.3730  0.0700 

0.0042 

27 

0.918 

-0.6334 

0.2597  0.0502 

0.0018 

29 

0.965 

-0.3072 

0.1191  0.0278 

0 . 0005 

30 

SUM 

0.988 

-0.0980 

-0.1762 

0.0370  0.0069 

0.0842  0.1157 

0.0009 

BKT, VBKT, 

HKT, VHKT 

0.1332  -0.0065 

0.0015 

BKQ, VBKQ, 

HKQ , VHKQ 

0.0310  0.0039 

0.0000 

DTMB 

DRAG  KT 

0.1333  KQ-  0.03203 

ETAO- 

Cf«0, 

.0000  KT 

-  0.1347  KQ«  0.03189 

ETAO- 

o 

a 

u 

.0035  KT 

«  0.1262  KQ-  0.03597 

'.  ETAO- 

Cf-0. 

.0070  KT 

-  0.1177  KQ*  0.04005 

ETAO- 

Cf-0 . 

.0100  KT 

-  0.1104  KQ-  0.04355 

ETAO- 

CDrag 

vx 

VR 

VT 

UANW  UTNW 

UAUW  UTUW 

0.0070 

1.050 

0.000 

-0.209 

-0.050  -0.026 

0.139  -0.273 

0.0000 

1.050 

0.000 

-0.186 

-0.039  -0.076 

0.105  -0.162 

0.0000 

1.050 

0.000 

-0.151 

0.037  -0.172 

0.143  -0.110 

0.0000 

1.051 

0.000 

-0.131 

0.140  -0.239 

0.229  -0.165 

0.0000 

1.053 

0.000 

-0.128 

0.214  -0.252 

0.248  -0.184 

0.0000 

1.056 

0.000 

-0.129 

0.230  -0.210 

0.178  -0.125 

0.0000 

1.058 

0.000 

-0.126 

0.179  -0.137 

0.092  -0.052 

0.0000 

1.058 

0.000 

-0.125 

0.127  -0.097 

0.068  -0.032 

0.0000 

1.058 

0.000 

-0.123 

0.057  -0.060 

0.069  -0.033 

0.0000 

1.057 

0.000 

-0.120 

-0.032  -0.028 

0.106  -0.066 

(SUCTION 

-  0) 

(PARTIAL 

TABLE) 

100 

100 

BKMZ 

CIRC 

0.0015 

-0.2212 

0.0012 

-0.0087 

0.0301 

0.4086 

0.0898 

0.8824 

0.1761 

1.3206 

0.2819 

1.6760 

0.3985 

1.9292 

0.5160 

2.0725 

0.6230 

2.1076 

0.7062 

2.0293 

0.7518 

1.8547 

0.7447 

1.5878 

0.6698 

1.2355 

0.5172 

0.8120 

0.2770 

0.2937 

0.1096 

-0.0398 

-0.0020 

0.0002 

0.7288 

0.7394 

0.6142 

0.5144 

0.4438 


* 


♦ 


1 
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